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ABSTRACT 

Particle acceleration at plasma shocks appears to be ubiquitous in the uni- 
verse, spanning systems in the heliosphere, supernova remnants, and relativistic 
jets in distant active galaxies and gamma-ray bursts. This review addresses some 
of the key issues for shock acceleration theory that require resolution in order to 
propel our understanding of particle energization in astrophysical environments. 
These include magnetic field amplification in shock ramps, the non-linear hydro- 
dynamic interplay between thermal ions and their extremely energetic counter- 
parts possessing ultrarelativistic energies, and the ability to inject and accelerate 
electrons in both non-relativistic and relativistic shocks. Recent observational 
developments that impact these issues are summarized. While these topics are 
currently being probed by astrophysicists using numerical simulations, they are 
also ripe for investigation in laboratory experiments, which potentially can pro- 
vide valuable insights into the physics of cosmic shocks. 

1. Introduction 



Supersonic flows abound in the cosmos, as do emission regions exhibiting non-thermal 
radiation. The intimate connection between the two establishes that particle acceleration in 
astrophysical shocks is germane to many systems, ranging from the heliosphere, to stars of 
various sorts expelling winds, to supernova remnants, to extragalactic jets and gamma-ray 
bursts. In the case of heliospheric shocks such as travelling interplanetary discontinuities and 
planetary bow shock environs, we can immerse ourselves in the plasma experiment via in 
situ spacecraft measurements of non-thermal ions, electrons and turbulent magnetic fields. 
While localized and therefore sparse in terms of the spatial sampling, these observations 
do provide profound insights into the complexity of the shock acceleration phenomenon. In 
astrophysical sites beyond the solar system, our role is passive, as observers of signals from 
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remote sites of acceleration. Moreover, the information on plasma properties is subject to 
a convolution with radiative processes, complicated source morphology within our spatial 
resolution scale, and propagational modification along the line of sight to sources. 

Observationally, radio, optical and X-ray telescopes have provided groundbreaking in- 
sights into the shock acceleration phenomenon, due to advances in the angular resolution 
and spectral sensitivity. In addition, the gamma-ray field is generating a greater understand- 
ing of source energetics as we find that many non-thermal astronomical sources emit most 
of their power in the gamma-rays. Progress on the theoretical front builds on the obser- 
vational advances, and has turned more to computer simulations due to the dramatically 
enhanced speed of computers over the last two decades. Yet astrophysical code verification 
is an increasingly salient issue as their complexity and computational demands burgeon. 
This provides a niche for laboratory plasma experiments that are tailored for the problem 
of astrophysical particle acceleration. Supersonic flows can be generated in controlled envi- 
ronments, and work along these lines focuses naturally on using lasers to mimic blast waves 
and jets, supernovae and supernova remnants (e.g. Borovsky et al. 1984; Drake et al. 1998; 
Shigemori, et al. 2000; Kang et al. 2001; Lebedev et al. 2002; Woolsey, Courtois & Dendy 
2004), and probe hydrodynamic and magnetohydronamic aspects. The question of scalabil- 
ity of such laboratory findings to astrophysical systems (e.g. Ryutov et al. 2001; see also 
the review of Remington, Drake & Ryutov 2006) is obviously of central importance. Here, 
an offering on some topical issues for particle acceleration at astrophysical shocks is made, 
to provide a basis for the community in high energy density plasma physics to help identify 
germane astrophysical problems that might be well suited for interdisciplinary investigation. 



2. Magnetic Field Enhancements in Shocks 

One of the key properties of shock structure that is germane to the acceleration of high 
energy cosmic rays is the strength of the magnetic field B near the shock. In astrophysi- 
cal shocks this cannot really be measured directly, since there is generally a lack of viable 
spectral line diagnostics: the Zeeman effect and cyclotron emission/absorption features are 
generally broadened, small or non-existent in diffuse, turbulent shock environs. Normally, 
proximity of an emission region to a stellar surface, such as in white dwarfs and neutron 
stars is required to afford precise magnetic field measurements. In heliospheric shocks, mag- 
netometer data discern the chaotic nature of pre- and post-shock fields, and a prominent 
property appears to be (e.g. Baring et al. 1997) a shock-induced compression of the field in 
the downstream region that is fairly close to magnetohydrodynamic (MHD) determinations 
that are derived from momentum and energy flux conservation across the shock, i.e. the 
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so-called Rankine-Hugoniot conditions (e.g. Drury 1983; Jones & Ellison 1991). If the in- 
teraction of charged particles with shock-associated field turbulence is gyroresonant at the 
Doppler-shifted cyclotron frequency (e.g., see Melrose 1980), as is expected for Alfven and 
whistler modes, then the acceleration timescale r naturally scales as the gyroperiod (i.e., 
t ~ oc 1/B; see Forman, Jokipii & Owens 1974; Drury 1983) and the corresponding 
diffusive lengthscale is comparable to the Larmor radius r g ( oc 1/B ). Hence the magnitude 
B establishes the temporal, spatial and energy scales of acceleration at a shock, and so is a 
critical parameter for the energization process. 

Since the general paradigm of galactic cosmic rays (CRs) invokes supernova remnants 
(SNRs) as the sites for their production (see Drury 1983 for a review), knowledge of the 
field strength in proximity of their shocks is vital. Directional information can be obtained 
on fairly large spatial scales via radio synchrotron polarization data (e.g. see Rosenberg 
1970; Downs & Thompson 1972; Anderson, Keohane & Rudnick 1995, for Cassiopeia A), 
but values of B = |B| are not forthcoming. Estimates for the field strength can be inferred 
by modeling the continuum flux level in a given wavelength band, but these are subject 
to a number of assumptions about the medium, for example the mean density n e and the 
Lorentz factors 7 e of the radiating electrons. The understanding of the character of shells 
and interiors of SNRs has recently been advanced by groundbreaking observations with the 
Chandra X-ray Observatory, enabled by its impressive angular resolution coupled with its 
spectral capabilities. Of particular interest is the observation of extremely narrow non- 
thermal (1.2 - 2.0 keV) X-ray spatial profiles in selected remnants (see Long et al. 2003 
and Bamba et al. 2003 for the northeast limb of SN1006; Vink & Laming 2003 for Cas A; 
for theoretical modeling see Ellison & Cassam-Chenai 2005, and Volk et al. 2005), typically 
less than 5-10 arcsec across. Upstream of these shell shocks, the X-ray emission, which 
is thought to be synchrotron in origin, drops to effectively zero. These strong brightness 
contrasts between the shell, and the outer, upstream zones correspond to flux ratios exceeding 
R ^ 50 . The narrowness of profiles along image scans argues for the shocks being aligned 
perpendicular to the sky, i.e., offering no projectional smearing in the images. Note also 
that the surface brightness angular profiles in SN1006 and Cas A are much broader for the 
thermal X-rays (0.5-0.8 keV) and the radio synchrotron than for the non-thermal X-rays. 

If the synchrotron mechanism is indeed responsible for non-thermal Chandra emission, 
the electrons contributing to the Chandra signal are probably in a strongly-cooling regime: 
see Baring et al. (1999) for a comprehensive discussion of SNR cooling parameter space. 
Since the synchrotron cooling rate for an electron scales as ^B 2 , then the flux ratio R is 
approximately a measure of the ratio of B 2 downstream (d) to upstream (u). The observed 
lower bounds to R considerably exceed values R ^ 16 expected for magnetohydrodynamic 
compression at the shocked shell; at a plane-parallel shock with B along the shock normal 
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there is no field compression, while in a strong (i.e. high sonic Mach number) perpendic- 
ular shock with B in the shock plane, B^jB u ~ 4. Hence, the pronounced brightness 
contrast is taken as strong evidence of magnetic field amplification in the shock pre- 
cursor/ramp upstream. Higher fields are obviously advantageous to cosmic ray production 
issues. Historically-accepted values of B ~ 1-10/x Gauss (i.e. 0.1-1 nanoTesla) are some- 
what too small to permit acceleration in SNR shocks of ages around 10 3 - 10 4 years right up 
to the cosmic ray knee at ~ 3 x 10 15 eV (e.g. see Lagage & Cesarsky 1983). This problem 
has spawned the suggestion (Jokipii 1987) that relatively ineffective diffusive transport of 
particles orthogonal to the mean field direction in quasi-perpendicular regions of SNR shocks 
can speed up acceleration of ions to higher energies, helping access the knee. Yet, this en- 
hanced rapidity is accompanied by reduced efficiency of cosmic ray injection from thermal 
energies (Ellison, Baring & Jones 1995). Hence, truly larger fields provide a cleaner path for 
acceleration in remnants to reach the cosmic ray knee (e.g. Kirk & Dendy 2001). 

The X-ray observational developments have been accompanied by theoretical proposals 
of magnetic field amplification in the upstream shock precursor. Most notable has been 
the work of Lucek & Bell (2000), and subsequent papers such as Bell (2004), Amati & 
Blasi (2006), and Vladimirov, Ellison & Bykov (2006). The idea of Lucek & Bell is that high 
energy cosmic rays in strong shocks could amplify B when streaming upstream, adiabatically 
transferring energy to the turbulent field by pushing against it, simultaneously decelerating 
the upstream flow. If this process is efficient, the rate of work done on the upstream Alfven 
turbulence of energy density U A naturally scales roughly with the CR pressure gradient: 
dU A /dt = WaIVPcrI • Here v A = B/y/Atrp is the Alfven speed, and Pgr is the cosmic ray 
pressure. The associated field amplification should then scale as (5B/B) 2 ~ M a Pcr/ pu 2 in 
an upstream flow of speed u u and mass density p ; this then becomes very effective for high 
Alfvenic Mach number (i.e. M. A = u u /v A ^> 1), strong shocks that generate large cosmic 
ray pressures. While this hypothesis is reasonable, demonstrating it is non-trivial. Various 
MHD-type simulations have been employed by Bell, such as in Bell (2004), where large-scale 
currents are used to drive instabilities that amplify the upstream field. The persistence 
of currents on large scales is unclear, particularly due to the action of Debye screening. 
Moreover, self-consistent physical connection between the cosmic rays of large Larmor radii 
and the field turbulence of much shorter wavelengths is extremely difficult to explore with 
MHD or plasma simulations, due to the wide disparity in spatial scales involved. This is 
an issue also for a growing number of particle-in-cell (PIC) simulations (e.g. Silva et al. 
2003; Hededal et al. 2004; Nishikawa et al. 2005; see Section @] below) used to explore field 
enhancement via the Weibel instability in relativistic shocks; such developments are not that 
salient for the problem of amplifying Alfven turbulence in non-relativistic shocks, and mostly 
probe the inertial scales of thermal ions and electrons defined by their plasma frequencies. 
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3. Non-linear Feedback Between the Acceleration and the Hydrodynamics 

Non-relativistic collisionless shocks can be highly efficient accelerators, placing 10-50% 
of the bulk flow kinetic energy into non-thermal particles. Evidence from theory, computer 
simulations, and spacecraft observations supports this conclusion; in particular, see Ellison, 
Mobius, & Paschmann (1990), for a study of the Earth's bow shock, and Drury (1983), 
Blandford & Eichler (1987), and Jones & Ellison (1991) for reviews. With such efficiencies, 
the accelerated particles acquire a sizable fraction of the total energy budget, influencing 
the shock hydrodynamics, and therefore also the fraction of energy going into accelerated 
particles, in a non-linear manner. The modified flow velocity spatial profile in the shock 
deviates from the familiar step-function form in test-particle acceleration scenarios, with the 
energetic particles pushing against the upstream flow and decelerating it far ahead of the 
shock discontinuity. Accordingly an upstream shock precursor forms, with declining flow 
velocity as the shock is approached. This structure alters the shape of the energetic particle 
distribution from a power-law in momentum (e.g., Eichler 1984; Ellison & Eichler 1984; 
Ellison, Baring & Jones 1996; Berezhko et al. 1996; Malkov 1997; Blasi 2002), the canonical 
test-particle form where the diffusively-accelerated particle distribution samples no spatial 
or momentum scale. The index a — (r + 2)/(r — 1) of this power-law dn/dp oc p~° is 
purely a function of the compression ratio r = u u /ud of upstream ( u u ) to downstream ( w d ) 
flow speed components normal to the shock in the shock rest frame (e.g. see Drury 1983; 
Jones & Ellison 1991), and is independent of the magnetic field orientation or the nature 
and magnitude of the turbulence effecting diffusive transport in the shock neighborhood. 

The spatial variation of the upstream flow in strong shocks that are efficient accelerators 
eliminates the scale independence. Since the highest energy particles have greater diffusive 
mean free paths A (generally true for gyroresonant interactions with MHD turbulence, and 
certainly so near the Bohm diffusion limit A ~ r g ), they diffuse farther into the upstream 
shock precursor against the convective power of the flow, and therefore sample greater effec- 
tive velocity compression ratios r . Accordingly, they have a flatter distribution, yielding a 
distinctive concavity to the overall particle spectrum, i.e. a is now a declining function of 
momentum p . These departures from power-law behavior amplify the energy placed in the 
particles with the greatest momenta, which in turn feeds back into the shock hydrodynamics 
that modify the spatial flow velocity profile. Traveling discontinuities possessing this complex 
feedback are termed non-linear shocks, the non-linear label being ascribed to the inter- 
play between the macroscopic dynamics and the microscopic acceleration process. Clearly, 
the possible magnetic field amplification in the upstream precursor that was discussed in 
Section [2] contributes to the overall dynamics/energy budget of the magnetohydrodynamic 
flow, and so intimately influences this non-linear aspect of astrophysical shocks. 
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The deviations from power-law distributions obviously impact the radiation signatures 
produced by these particles, with alterations in the fluxes expected in X-ray and TeV gamma- 
ray bands in remnants, differing by as much as factors of 3-10 from traditional test-particle 
predictions (e.g., see Baring, et al. 1999; Ellison, Slane & Gaensler 2001; Berezhko et al. 
2002; Baring, Ellison & Slane 2005). Conclusively confirming the existence of this non-linear 
spectral concavity is a major goal that is inherently difficult, since it demands broad, multi- 
wavelength spectral coverage. There is a limited suggestion of concavity in radio data for 
Tycho's and Kepler's SNRs (Reynolds & Ellison 1992), and in a multi- wavelength modeling 
of SN 1006 (Allen, Houck & Sturner 2004; see also Jones et al. 2003 for inferences from 
radio and infra-red data from Cassiopeia A), but this task really looks ahead to the launch 
of the GLAST gamma-ray mission in late 2007, when, in conjunction with ground-based 
Atmospheric Cerenkov Telescopes probing the TeV band, it may prove possible to determine 
gamma-ray spectra from SNRs spanning over 3 decades in energy. 

In the meantime, an interesting astrophysical manifestation of these non-linear effects 
has been offered by SNR observations by the Chandra X-ray Observatory, looking instead 
at the thermal populations. Inferences of ion temperatures in remnant shocks can be made 
using proper motion studies, or more direct spectroscopic methods (e.g. Ghavamian et al. 
2003). For the remnant IE 0101.2-7129, Hughes et al. (2000) used a combination of ROSAT 
and Chandra data spanning a decade to deduce an expansion speed. Electron temperatures 
T e are determined by line diagnostics, via both the widths, and the relative strengths for 
different ionized species. From these two ingredients, Hughes et al. (2000) observed that, 
in selected portions of the SNR shell, 3£;T e /2 < 3£;T p /2 ~ m p (3w u /4) 2 /2 . Therefore, the 
electrons were considerably cooler than would correspond to equipartition with thermal 
protons heated in a strong shock with an upstream flow speed of u u : the thermal heating 
is assumed comparable to the kinematic velocity differential u u — ~ 3u u /4. The same 
inference was made by Decourchelle et al. (2000) for Kepler's remnant, and by Hwang et al. 
(2002) for Tycho's SNR. This property of comparatively cooler electrons may be indicative of 
them radiating very efficiently radiating. Or it may suggest that the protons are cooler (i.e. 
3kT p /2 <C m p (3u u /A) 2 /2 ) than is widely assumed in the test-particle theory, the conclusion 
drawn by Hughes et al. (2000) and Decourchelle et al. (2000). This effect is naturally 
expected in the non-linear shock acceleration scenario: as the highest energy particles tap 
significant fractions of the total available energy, they force a reduction in the thermal gas 
temperatures. Such feedback can profoundly influence shock layer thermalization, inducing 
significant interplay with electrostatic equilibration between low energy electrons and ions, 
an issue addressed in these proceedings by Baring & Summerlin (2006). Note that non-linear 
modifications may vary strongly around the shocked shell of an SNR, since the obliquity angle 
© Bn ,u of B to the shock normal varies considerably between different rim locales. 
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4. The Character of Relativistic Shocks 

Relativistic shocks, for which the upstream flow Lorentz factor 7 n = 1/ a/1 — (u u /c) 2 
considerably exceeds unity, are less well researched than their non-relativistic counterparts, 
not in small part due to their greater cosmic remoteness: they predominantly arise in ex- 
tragalactic locales like jets in active galaxies, and gamma-ray bursts. Yet, because of such 
associations, they are now quite topical. Diffusive test-particle acceleration theory in par- 
allel, (i.e., 6 B „, U = 0°) relativistic shocks identifies two notable properties in such systems: 
(i) particles receive a large energy kick AE ~ 7 u ^c 2 in their first shock crossing (e.g., 
Vietri 1995), but receive much smaller energy boosts for subsequent crossing cycles (factors 
of around two: e.g., Gallant & Achterberg 1999; Baring 1999); (ii) a so-called 'universal' 
spectral index, a ~ 2.23 exists in the two limits of 7 U ^> 1 and small angle scattering, i.e., 
59 < 1/7 U (e.g., Kirk et al. 2000; see also Bednarz & Ostrowski 1998; Baring 1999; Ellison 
and Double 2004). Here, 59 is the average angle a particle's momentum vector deviates in 
a scattering event, i.e. an interaction with magnetic turbulence. 

These characteristics are modified in parallel, mildly relativistic shocks with 7„ ~ 1 . 
In such shocks, the distribution dn/dp remains a power-law (scale- independence persists), 
but hardens (a decreases) as either 7„ drops, or the scattering angle, 59, increases (e.g., 
Ellison, Jones & Reynolds 1990; Baring 1999; Ellison & Double 2004; Baring 2004), even if 
the compression ration r = u u /ud is held constant (it usually increases with declining 7 U 
due to a hardening of the Juttner-Synge equation of state). These effects are consequences 
of large kinematic energy kicks particles receive when scattered in the upstream region after 
transits from downstream of the shock. It is particularly interesting that when scattering 
conditions deviate from fine pitch-angle-scattering regimes with 59 <C l/7„ , the power-law 
index is dependent on 59 , with a continuum of spectral indices being possible (Ellison & 
Double 2004; Baring 2004). Then the nature of the turbulence is extremely influential on 
the acceleration outcome, so that understanding the turbulence is of paramount importance. 
This sensitivity of a to the field fluctuations when 59 ^ l/7«, a large angle scattering 
domain, contrasts the canonical nature of a in non-relativistic shocks mentioned above. 

In jets and gamma-ray bursts, ultra- relativistic shocks are typically highly oblique due 
to the Lorentz transformation of ambient, upstream magnetic fields to the shock rest frame. 
This introduces an added dimension of variation, with increasing Bn , u dramatically steep- 
ening the power-law, i.e. increasing a . This is naturally expected since such systems are 
highly superluminal, that is, there exists no de Hoffman- Teller (1950) shock rest frame where 
the flow velocities are everywhere parallel to the mean magnetic field (which would corre- 
spond to large scale electric fields being zero everywhere). Therefore, relativistic shocks 
are much less efficient accelerators because particles convect more rapidly away downstream 
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from the shock (e.g. Begelman & Kirk 1990). In oblique, relativistic shocks, a , and indeed 
the efficiency of injection from the thermal particle population, also depend on the ability 
of turbulence to transport particles perpendicular to the mean downstream field direction 
(Ellison & Double 2004; Niemiec & Ostrowski 2004). This perpendicular transport cou- 
ples directly to the magnitude (5B/B) 2 and power spectrum of field fluctuations, i.e. the 
strength of the scattering. Steep spectra (<7 ^ 4) result unless the ratio of the diffusive 
mean free paths perpendicular to and parallel to B is comparable to unity, which defines the 
Bohm diffusion regime. In summation, for relativistic shocks, the spectral index is sensitive 
to the obliquity 6 B n, u of the shock, the nature of the scattering, and the strength of the 
turbulence or anisotropy of the diffusion. These properties are reviewed in Baring (2004). 

Observational vindication of these theoretical predictions is clearly mandated. This 
is not readily forthcoming, since the only accessible information involves a convolution of 
shock acceleration and radiation physics. Yet, it is clear, for example in gamma-ray bursts 
(GRBs), that data taken from the EGRET experiment on the Compton Gamma- Ray Ob- 
servatory (CGRO) suggest a broad range of spectral indices (Dingus 1995) for the half dozen 
or so bursts seen at high energies. This population characteristic is commensurate with the 
expected non-universality of o just discussed. Yet it is important to emphasize that the 
power-law index is not the only acceleration characteristic germane to the GRB problem: 
the shapes of the particle distributions at thermal and slightly suprathermal energies are 
also pertinent. This energy domain samples particle injection or dissipational heating in the 
shock layer, and is readily probed for electrons by the spectrum of prompt GRB emission by 
the BATSE instrument on CGRO. Tavani (1996) obtained impressive spectral fits to several 
bright BATSE bursts using a phenomenological electron distribution and the synchrotron 
emission mechanism. While there are issues with fitting low energy (i.e. ^ 100 keV) spectra 
in about 1/3 of bursts (e.g. Preece et al. 1998) in the synchrotron model, this radiative 
mechanism still remains the most popular candidate today for prompt burst signals. 

Tavani's work was extended recently by Baring & Braby (2004), who provided addi- 
tional perspectives, using acceleration theory to underpin a program of spectral fitting of 
GRB emission using a sum of thermal and non-thermal electron populations. These fits 
demanded that the preponderance of electrons that are responsible for the prompt emis- 
sion constitute an intrinsically non-thermal population. That is, the contribution to 
the overall electron distribution that comes from a Maxwell-Boltzmann distribution is com- 
pletely dominated by a non-thermal contribution that, to first order, can be approximated 
by a power-law in energy truncated at some minimum electron Lorentz factor. This re- 
quirement of non-thermal dominance strongly contrasts particle distributions obtained from 
acceleration simulations, as is evident in a host of the references cited on acceleration the- 
ory above: the non-thermal particles are drawn directly from a thermal gas, a virtually 
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ubiquitous phenomenon. This conflict poses a problem for acceleration scenarios unless (i) 
radiative efficiencies for electrons in GRBs only become significant at highly superthermal 
energies, or (ii) shock layer dissipation in relativistic systems can suppress thermalization of 
electrons. A potential resolution to this dilemma along the lines of option (i) is that strong 
radiative self-absorption could be acting, in which case the BATSE spectral probe is not 
actually sampling the thermal electrons. It is also possible that other radiation mechanisms 
such as Compton scattering, pitch- angle synchrotron, or jitter radiation may prove more 
germane. Discerning the radiation mechanism(s) operating in bursts is a foremost goal of 
future research, and will be facilitated by the GLAST mission, with its good sensitivity in 
the 5 keV - 300 GeV band, in conjunction with NASA's current GRB flagship venture, Swift. 

Option (ii) is a conjecture that has no definitive simulational evidence to support it 
at present. The most comprehensive way to study dissipation and wave generation in col- 
lisionless shocks is with PIC simulations, where particle motion and field fluctuations are 
obtained as solutions of the Newton-Lorentz and Maxwell's equations. Rich in their tur- 
bulence information, these have been used extensively in non-relativistic, heliospheric shock 
applications, and more recently, relativistic PIC codes have blossomed to model shocks in 
various astrophysical systems. PIC simulation research has largely, but not exclusively, fo- 
cused on perpendicular shocks, first with Gallant et al.(1992), Hoshino et al. (1992), and 
then Smolsky & Usov (1996), Shimada & Hoshino (2000), Silva et al. (2003), Nishikawa 
et al. (2003, 2005), Spitkovsky & Arons (2004), Hededal et al. (2004), Liang & Nishimura 
(2004), Medvedev et al. (2005) and Hededal & Nishikawa (2005). These works have explored 
pair shocks, ion-doped shocks, Poynting flux-dominated outflows, and low-field systems with 
dissipation driven by the Weibel instability, in applications such as GRBs and pulsar wind 
termination shocks. PIC simulations are dynamic in nature, and rarely achieve a time- 
asymptotic state. Even in the minority of cases where there is some evidence of acceleration 
beyond true thermalization, none of these works has demonstrated the establishment of an 
extended power-law that is required in modeling emission from GRBs and active galactic 
nuclei. This is perhaps due to the severely restricted spatial and temporal scales of the simu- 
lations, imposed by their intensive CPU and memory requirements; these limit the modeling 
of realistic electron-to-proton mass ratios, full exploration of three-dimensional shock physics 
such as diffusive transport, and addressing the wide range of particle momenta encountered 
in the shock acceleration process. In particular, it is difficult to establish a broad inertial 
range for cascading MHD turbulence when the maximum spatial scales in the simulation are 
not orders of magnitude larger than the principal ion inertial scales. Definitively observing 
the injection of electrons from a thermal population to establish a truly non-thermal distri- 
bution remains a pressing goal of plasma simulations, both for relativistic shocks and, as has 
been understood for more than two decades, also their non-relativistic cousins. 
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5. Conclusion 

This review is by no means a complete presentation of the topical issues for the shock ac- 
celeration problem, but it does offer a fair sampling suitable for motivating interdisciplinary 
activity. It is clear that several issues could benefit substantially from input from laboratory 
experimentation on the high energy density physics/astrophysics interface. One key question 
is whether or not ambient magnetic fields are amplified by both non-relativistic and rela- 
tivistic shocks beyond standard MHD expectations. If so, is the amplification electrostatic in 
origin, or is it connected to energetic particles accelerated by the shock? It would be impor- 
tant to discern whether there are differences between high and low Alfvenic Mach number 
systems, i.e. what role the ambient magnetic field plays in controlling the outcome. Another 
question concerns whether or not suprathermal electrons and ions can actually be seen, and 
whether one can identify their origin. It is salient to ascertain if they are diffuse in nature, 
or if they form coherent beams, both of which are seen at traveling shocks embedded in the 
solar wind. Also, if acceleration is observed, then identifying the role the highest energy 
particles have in modifying the shock hydrodynamics and the thermal structure of the shock 
layer would help solve an outstanding problem that has long been a principal goal within 
the cosmic ray community. Finally, specifically concerning relativistic systems, it would be 
desirable to elucidate how the distributions of any accelerated particles seen depend on ex- 
ternal quantities such as the field obliquity and speed of the shock, whether thermal electrons 
can be suppressed relative to accelerated ones, and if there is an identifiable connection with 
the field turbulence near the shock. These are demanding goals, yet terrestrial experiments 
are very useful for probing global aspects of shock problems, and in particular for extracting 
insights into hydrodynamic and MHD behavior. In order to make progress, it is essential 
to prepare an experimental setup that is as tenuous as possible, to mimic the collisionless 
(in the Spitzer sense) shock environments offered throughout the cosmos. At this juncture, 
exciting prospects are on the horizon for this interdisciplinary forum, with contributions to 
be found in the laboratory, in computer simulations, and in astronomical observations, all 
of which can benefit from cross-fertilization with each other. 
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